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Numerical calculations of the average dipole-coupling energy Edi p in two-dimensional disordered 
magnetic nanostructures are performed as function of the particle coverage C. We observe that 
-Edip scales as Edi p oc C a with an unusually small exponent a* ~ 0.8-1.0 for coverages C < 20%. 
This behavior is shown to be primarly given by the contributions of particle pairs at short distances, 
which is intrinsically related to the presence of an appreciable degree of disorder. The value of a* 
is found to be sensitive to the magnetic arrangement within the nanostructure and to the degree of 
disorder. For large coverages C > 20% we obtain -Edip oc C a with a = 3/2, in agreement with the 
straighforward scaling of the dipole coupling as in a periodic particle setup. Taking into account the 
effect of single-particle anisotropies, we show that the scaling exponent can be used as a criterion 
to distinguish between weakly interacting (a* ~ 1.0) and strongly interacting (a* ~ 0.8) particle 
ensembles as function of coverage. 

PACS numbers: 75.75.+a, 75.70.Ak, 75.50.Lk, 61.46.+W 



I. INTRODUCTION 



Nanostructured materials consisting of interacting 
magnetic particles are currently the subject of an intense 
research activity, driven by their fundamental interest 
and technological perspectives. Numerous experimental 
and theoretical studies have been performed for various 
two-dimensional (2D) 1 and three-dimensional (3D) 2 - 3 - 4 - 5 
arrangements having different degrees of structural and 
magnetic disorder. The magnetic behavior of these sys- 
tems is determined both by single-particle properties, like 
lattice or shape anisotropies, and by interparticle interac- 
tions, like magnetic dipole coupling or indirect exchange. 
The relative importance of these two types of contribu- 
tions can be tuned experimentally, at least to some ex- 
tent, by varying the particle-size distribution, the nanos- 
tructure morphology, and the average interparticle dis- 
tance, for example. For the weakly interacting regime 
the single-particle properties dominate, and the inter- 
actions can be handled as perturbations. However, in 
the most interesting case of strongly interacting systems, 
where the interparticle couplings are stronger than the 
single-particle energies, an explicit treatment of the in- 
teractions is crucial. Taking into account the complexity 
of the problem, particularly in the presence of disorder, it 
is a challenge to derive simple trends for the dependence 
of the interaction energy on the main nanostructure char- 
acteristics. 

In this study we focus on the dipole coupling energy 
.Edip for magnetic particle ensembles with varying par- 
ticle coverages C and different kinds of nanostructures^ 
Since the dipole interaction energy between two particles 
scales with the average interparticle distance R as R~ 3 , 
and since C is propotional to R~ D , D being the spatial 
dimension, one could in principle expect that the average 
dipole energy -Edip (C) should be proportional to C a , with 
a = 3/-D. However, several previous studies indicate a 



different scaling behavior for disordered dipole-coupled 
systems i&i 5 . For example, a weaker decrease of -Edip char- 
acterized by an effective exponent a* — 0.8 ± 0.1 has 
been found experimentally upon thinning of ferrofluids 
(D = 3). 3 In addition, the strength of the dipole coupling- 
seems to be underestimated as compared to measure- 
ments, when one calculates its effect on the energy barrier 
for single-particle magnetic reversals^ Moreover, numer- 
ical results for -Edip of a ferrofluid monolayer (D = 2) 
with different particle coverages reveal significant devi- 
ations from the a = 3/2 scaling* 5 . These experimental 
and theoretical findings suggest that the effective dipole 
energy in disordered systems decreases more slowly with 
decreasing coverage than originally expected. 

The purpose of this study is to determine the role 
of global nanostructure characteristics, such as coverage 
and degree of disorder, on the average dipole-coupling 
energy of 2D particle ensembles. In particular we show 
that the scaling exponent a, -Ed; p oc C Q , allows to dis- 
cern whether £di P is dominated by the average inter- 
particle distance R, as in a periodic particle setup, or 
by short-distance contributions. Moreover, taking into 
account single-particle anisotropies we demonstrate that 
the value of a can be used as a criterion to distinguish 
between weakly and strongly interacting disordered par- 
ticle ensembles. 



II. MODEL SIMULATIONS 

The 2D nanostructure is modeled using a square unit 
cell containing non-overlapping, sphere-shaped magnetic 
particles, and periodic boundary conditions. We consider 
disordered particle ensembles with a fully random distri- 
bution within the unit cell, as well as disturbed square 
arrays where the particles are scattered around periodic 
lattice sites by applying a Gaussian distribution P(r) 
with standard deviation <tr. For such planar particle 



ensembles the overall coverage C is related to the aver- 
age interparticle distance R > 2r by C = it (r Q /R) 2 . A 
uniform-particle size distribution with N ~ 24000 atoms 
per particle is assumed, which corresponds to a particle 
radius ro ~ 16 do, where ao is the nearest-neighbor inter- 
atomic distance. Due to the strong interatomic exchange 
interaction and for particle diameters smaller than typ- 
ical magnetic domain wall widths, each particle can be 
regarded as a single-domain Stoner-Wohlfarth magnet- 
carrying a giant spin M ~ N fi at , where /i at is the mag- 
netic moment per atom. For simplicity, the particle mag- 
netizations Mj are restricted to be within the plane of 
the nanostructure. The magnetic configuration is thus 
characterized by the set of in-plane angles {4>i}, with 
i = 1,...,N. For each given particle distribution and 
magnetic configuration we determine the magnetic en- 
ergy 

E = E anis + E dip = -K 2 J2( M i n *) 2 

i 

+\ [M( M j - 3(r« M,)(r tJ M.) r^] . (1) 

The first term represents the single-particle anisotropy 
energy of particle i, resulting from spin-orbit interac- 
tions and the shape anisotropy. For simplicity, a uni- 
axial anisotropy is assumed with strength K 2 per atom 
and easy axis directed along a given unit vector n^. The 
second term of Eq. is the dipole-dipole interaction be- 
tween magnetic moments Mj, where r%j = |ry| = — ij| 
is the distance between the centers of particles i and j. 
The infinite range of the dipole interaction is taken into 
account by applying an Ewald-type summation over all 
periodically arranged unit cells of the extended planar 
system^ In addition, we also include the leading cor- 
rection resulting from the dipole-quadrupole interaction, 
which takes into account the effect of finite particle sizes^ 
In order to determine the average magnetic energy 
E = -Eanis + Edip we start from an arbitrary initial guess 
of the angles <pi of each particle, and relax the magnetic 
configuration to a nearby local minimum of the total en- 
ergy, using a conjugated gradient method. The resulting 
relaxed states depend in general on the initial configu- 
rations and show strong magnetic noncollinearities. The 
large diversity of metastable states which is characteristic 
of disordered particle arrangements^^ is taken into ac- 
count by averaging E over many different random initial 
guesses (typically ~ 100) for the same spatial arrange- 
ment of the particles. Notice that this procedure is not 
intended to search preferently for the global energy min- 
imum. For comparison, E is also calculated for the same 
particle setup but assuming a parallel alignment of all 
particle magnetizations (i.e., a ferromagnetic state with 
4>i = 4>fm) and averaging over the common magnetization 
angle </>f m . Note that for a disordered particle ensemble 
such a ferromagnetic state does not usually correspond 
to a local energy minimum. In addition, a second aver- 
age over the spatial configuration of the nanostructure 
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FIG. 1: Average dipole energy |.Edi p (C)| per particle of ran- 
dom 2D ensembles of magnetic particles as function of the par- 
ticle coverage C, neglecting single-particle anisotropies (K2 = 
0). Results are given for relaxed magnetic states (circles), 
and for ferromagnetic ensemble magnetizations (crosses). The 
straight lines correspond to Edi p (C) oc C a , with a = 1 
(dashed) and a = 3/2 (dotted). 



is necessary in order to minimize spurious effects of the 
finite size of the unit cell. For this purpose, an appro- 
priate number (typically ~ 200) of different realizations 
of the unit cell is considered, using the same global vari- 
ables which characterize the particle ensemble. In the 
following calculations we use parameters appropriate for 
Fe, i.e., [x a t — 2.2 (13 and ao = 2.5 A, yielding the dipole 
energy unit E di = Hat/ a o — 0-19 K. Concerning the 
strength K2 of the uniaxial anisotropy we consider val- 
ues appropriate for films and particles of 3d transition 
metals (K 2 = 0.1-1 K/atom)ia 



III. RESULTS AND DISCUSSION 

First we focus on strongly interacting 2D systems of 
magnetic particles, neglecting the single-particle contri- 
bution (K2 — 0). In Fig.l results are given for the ab- 
solute value of the average dipole energy per particle 
|-E-di P (C)| as function of particle coverage C for a fully 
random particle ensemble. -Edip (C) I is shown for the 
relaxed magnetic arrangements as well as for the ferro- 
magnetic ones. From the log- log-plot of \Edi v {C)\ one 
observes that the exponent a — 3/2, corresponding to 
the straightforward scaling of the dipole interaction for 
D = 2, is approximately valid only for large coverages 
C > 20%. For smaller C, in the investigated range 
0.1% < C < 20%, one finds a different dependence with 
clearly smaller effective exponents. For the ferromag- 
netic arrangement we obtain af m ~ 1.0 and, quite inter- 
estingly, for the relaxed solution an even smaller value 
a* ~ 0.8 — 0.9. The numerical accuracy of the calculated 
exponents a^ m and a* is estimated to be about ±0.02. 

These remarkable results can be qualitatively ex- 
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plained by the following simple analysis. Consider 
two non-overlapping particles with radius r in a D- 
dimensional space coupled by an interaction E n (r) oc 
r~ n . Let the first particle be located at r — 0, and 
let P(r) be the probability to find the second one be- 
tween 2ro and a maximal distance i? max . The upper 
bound -R m ax is related to the particle concentration by 
C cx (ro/Rmax) , whereas the lower bound 2ro re- 
sults from the hard core (excluded volume) of the non- 
overlapping particles. The average interaction energy E n 
of such a particle pair is proportional to the average 



(r- n ) = / drr D - l P{r)r- 

J2r 
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drr u ~ v P(r). 

For a random distribution corresponding to a constant 
P(r) one obtains 
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for n = D. In the limit of low coverages, 2 r$ <C R n 
oo, this yields 
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for n > D, and a similar limit for n = D. For n < D 
one finds oc C n ^ D , which coincides with straight- 

forward scaling. Notice that for n > D the average in- 
teraction energy E n (R max ) for i? max — ► oo is no longer 
proportional to (-Rmax) Instead, an effective exponent 
n* = D is obtained. Consequently, as function of con- 
centration one obtains E n (C) oc C a with a* = 1. The 
physical reason behind this behavior is that only the up- 
per bound i?max scales with C, since it is proportional 
to the average interparticle distance, but not the lower 
bound 2rQ. Hence, the average energy is dominated by 
the interactions between particle pairs at distances of the 
order of the particle diameter, which is proportional to 
(2 ro)~", multiplied by the probability to find such a par- 
ticle pair, which is proportional to C . In other words, as 
a consequence of the relatively strong distance depen- 
dence of the interaction for n > D, in a diluted random 
system the average energy E n (C) is dominated by the 
short-distance contributions, and not by the energy cor- 
responding to the average interparticle distance. Notice 
that if n > D, E n oc C is independent of the particu- 
lar kind of interaction characterized by the number n. 
Moreover, this behavior of a random particle ensemble is 
also obtained for an exponentially decreasing interaction 
E oc cxp(— Ar) for all D. 

In the context of the present study corresponding to 
D = 2 and n — 3, one concludes that a 2D random en- 
semble of magnetic nanoparticles exhibits a stronger than 



expected average dipole interaction, which scales with the 
particle coverage as Edi P (C) oc C a with a* ~ 1.0. In 
contrast, if one varies the interparticle distance in a pe- 
riodic array the smallest distance between particle pairs 
also scales with coverage, yielding therefore the straight- 
forward scaling behavior with a = 3/D. Thus, by analyz- 
ing the exponent a one can distinguish whether the av- 
erage dipole energy is dominated by short-distance con- 
tributions or by the energy corresponding to the average 
interparticle distance. The difference between a and a* 
is significant for ID and 2D systems, whereas for 3D sys- 
tems one expects a ~ a* ~ 1. 

Another interesting effect is the difference in the ex- 
ponents a* and a^ m obtained for the relaxed and ferro- 
magnetic arrangements. This seems to be an indication 
of the formation of particle chains with correlated 'head- 
to-tail' alignments of the magnetizations of closeby parti- 
cles. Such magnetic rearrangements are possible only by 
relaxing the angles 4>i. In fact, recent calculations have 
shown that |i?dip | increases with increasing positional dis- 
order as a result of magnetic relaxations. ? Moreover, re- 
duced values of a* = 0.8 ±0.1 have also been observed 
in experiments on frozen high-density ferrofluids^ii 

By closer inspection of Fig.l one observes that the ex- 
ponent a* is not exactly the same for all coverages, since 
a slight change of slope appears in the curve of the re- 
laxed solution. Moreover, differences between a* and 
af m are found in the coverage range where magnetic re- 
laxation driven by the dipole coupling is significant, i.e., 
in the strongly interacting regime. This can be demon- 
strated by considering in addition to the dipole coupling 
a single-particle uniaxial anisotropy Ki with randomly 
distributed directions of the easy axes. Unless the av- 
erage anisotropy energy |i? an i s | is so strong that it dom- 
inates over |i?dip| even for large coverages, a crossover 
from the strongly interacting regime (|-E a nis| < |£-dip|) 
to the weakly interacting regime (|-E a nis| > l-^dipl) takes 
place with decreasing C at the crossover coverage C(Ki). 
This is shown in Fig. 2 where, as in Fig.l, |i?dip| and 
|-Eanis| are given as function of C for different values 
of K2- We have determined the relaxed magnetic ar- 
rangements in presence of local anisotropies, and have 
derived the scaling exponents a* corresponding to the 
dipole contribution of the magnetic energy. First of all, 
one observes that the scaling behavior, as estimated by 
Eq.© and depicted in Fig.l, is also valid in presence of 
single-particle anisotropies. In the coverage range where 
|-Edip| is dominated by short-distance contributions, i.e., 
for C < 20%, and in the strongly interacting regime 
C > C where |i?di P | > -Eanisl, a reduced value of a* < 1 
is clearly observed. As C is further decreased, a* changes 
gradually to a* ~ 1.0 below the crossover coverage C. 
This can be understood by recalling that magnetic re- 
laxation driven by dipole coupling is significant only in 
the strongly interacting regime. In the opposite case of 
weak interactions (C < C) relaxation due to dipole cou- 
pling is negligible, thus resulting in a* ~ aj- m ~ 1.0. 




FIG. 2: Average dipole energy |Z?di p (C)| for the relaxed mag- 
netic states as function of the particle coverage C for sev- 
eral values of the single-particle uniaxial anisotropy K2 = 
0.1 K/atom (diamonds), K2 = 0.2 K/atom (squares), and 
K 2 = 0.5 K/atom (triangles). For K 2 = 0.2 K/atom, the 
average anisotropy energy |-E an i s (C)| is also shown (dotted 
line, squares). The vertical lines indicate the anisotropy- 
dependent crossover coverages C^Kz), for which -Bdip = -E an is. 
For the sake of comparison the dashed curves show |_Edi P (C)| 
for K2 = for both the relaxed states (circles) and ferromag- 
netic states (crosses) as in Fig.l. 



FIG. 3: Scaling exponents a* and aJm °f the average dipole 
energy |i?dip(C')| oc C a corresponding to the relaxed and fer- 
romagnetic arrangements as function of the positional stan- 
dard deviation ur in units of the average interparticle distance 
R. A disturbed square particle array is assumed, neglecting 
single-particle anisotropics (K2 = 0). The exponents corre- 
spond to the coverage range 0.1% < C < 20%. The dashed 
curve refers to a* as derived from Eq.©. The horizontal lines 
indicate a* and af m calculated for a fully random particle en- 
semble. 



Note that for the largest considered value of K% the 
crossover happens directly from a = 3/2 to a* ~ 1.0, 
since here CiKi) ~ 20%. Consequently, the determi- 
nation of the scaling exponent a* allows one to discern 
between strongly and weakly interacting ensembles. 

The average anisotropy energy |-E an is| is of course in- 
dependent of coverage C in the absence of dipole inter- 
actions and other magnetic couplings. A dependence on 
C appears only if the interactions force the particle mag- 
netizations Mi to deviate from the local easy axes rij, cf. 
Eq.Q). Moreover, the exponent aif m for a ferromagnetic 
state is not affected at all by single-particle anisotropies. 
Since for such an arrangement a magnetic relaxation is 
not effective, the single-particle energy yields a contri- 
bution to the total magnetic energy which is indepen- 
dent of the coverage. For anisotropies with aligned easy- 
axes directions we obtain the same behavior of a* as 
for randomly distributed ones, since the dipole-coupling- 
induced relaxation of the particle magnetizations is not 
much affected by the easy axes distribution. 

In order to study the differences between a* and «f m 
in further detail we have also considered ID systems of 
dipole-coupled particles. In this case the particle mag- 
netizations are aligned parallel to the chain axes in the 
most stable state. A numerical investigation of such ID 
systems with randomly assembled particles yields effec- 
tive exponents aj m ~ a* ~ 1.0. This is quite reasonable 
since in this case the relaxed and ferromagnetic solutions 
coincide. In the case of ID stripes with finite widths, 



we observe a smooth crossover from a r * ~ 1 to a* ~ 0.8 
with increasing stripe width in the investigated coverage 
range 0.1% < C < 20%. These findings confirm the idea 
that the reduction of a* with respect to in 2D and 
3D nanostructures is the result of magnetic relaxations. 

To quantify the role of disorder on the coverage de- 
pendence of -Edip we have studied the crossover between 
periodic and fully random particle setups by considering 
systems with intermediate degrees of positional disorder. 
Here the single-particle anisotropy is neglected for sim- 
plicity {K2 = 0). Two-dimensional particle ensembles 
have been constructed with a quasi-periodic arrangement 
in which the particles are scattered around the sites of a 
square lattice following a Gaussian-like distribution with 
standard deviation or. Thus, a periodic square lattice 
corresponds to or = 0, and the particle arrangement be- 
comes increasingly disordered with increasing or. For 
each or the average dipole energy E^n p (C) has been de- 
termined as function of particle coverage C . The corre- 
sponding effective exponents a* and af m for the relaxed 
and ferromagnetic states are presented in Fig. 3 as func- 
tion of ctr/R, where R is the average interparticle dis- 
tance. For small or, i.e., < ctr/R < 0.2, one observes 
that a* ~ Q!f m ~ 3/2, in agreement with the straight- 
forward scaling of the dipole coupling. As or increases 
beyond orJR = 0.2, the exponents decrease rapidly, ap- 
proaching the limit of a fully random particle ensemble. 
In fact, previous calculations have shown that i?dip(cfl) is 
close to i?dip of random particle arrangements already for 
cr/R = 0.4 — O.61S These numerical results are corrobo- 
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rated by the simple analytical estimate given by Eq.@, 
using for P(r) a Gaussian distribution with mean i? max /2 
and standard deviation <tr. The different scaling proper- 
ties of ordered and disordered nanostructures confirm the 
importance of the sample characteristics on the effective 
scaling exponents. 

Finally, several other effects have been investigated 
which could in principle modify the values of the expo- 
nents a* . We have considered (i) unit cells with different 
numbers of particles, (ii) particle-size dispersions with 
Gaussian or log- normal distributions, (iii) disk-shaped 
particles with the same particle radius r$ as the pre- 
viously considered sphere-shaped particles, and (iv) the 
role of dipole-quadrupole corrections to the point-dipole 
sums. While the magnitude of Edi P is obviously affected 
by some of these contributions, in all cases the effective 
exponents a* and aj m remain unchanged within numer- 
ical accuracy. These results confirm the validity of the 
previously discussed trends and reinforce the idea that 
a* can be used as a parameter to characterize some as- 
pects of the magnetic behavior of interacting magnetic 
nanostructures with varying degrees of disorder. 

IV. CONCLUSION 

The average dipole energy i?di P of two-dimensional dis- 
ordered magnetic nanostructures has been determined by 
using numerical simulations and simple analytical esti- 
mates. A coverage dependence of the form \Edi P (C)\ oc 
C a is obtained with nearly constant a* over a wide 
coverage range. For 0.1% < C < 20% the exponent 
a* ~ 0.8-1.0 is significantly smaller than the value 



a = n/D = 3/2 expected from a straightforward scal- 
ing of the dipole interaction. This scaling results from 
dominating short-distance contributions, and is intrinsi- 
cally related to the presence of significant disorder. In 
contrast, for C > 20% the behavior of Edi P (C) is de- 
termined mainly by the average interparticle distance, 
which depends on coverage as R oc C~ X I D , and therefore 
yields a — 3/2, as for a periodic particle setup. A smooth 
transition between these two scaling regimes is observed 
with increasing degree of disorder. 

Taking into account the effect of single-particle 
anisotropics, we have shown that the value of the scaling 
exponent a* allows to distinguish between strongly and 
weakly interacting particle ensembles, since a* is sensi- 
tive to the magnetic arrangement within the nanostruc- 
ture. Indeed, one observes a crossover from a* ~ 0.8, 
when magnetic relaxations driven by dipole couplings 
dominate, to a* ~ 1.0, when the directions of the particle 
moments are determined mainly by local anisotropies. It 
would be very interesting to extend these investigations 
to 3D nanostructures in view of the numerous experi- 
mental results that are available for such systems. In ad- 
dition, the scaling behavior should be analyzed as func- 
tion of temperature, since it is likely to be affected by 
thermally induced crossings of energy barriers between 
different magnetic arrangements. 
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